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Abstract 



We study how some recently proposed noncontextuality tests based on 
quantum interferometry are affected if the test particles propagate as 
open systems in presence of a gaussian stochastic background. We show 
that physical consistency requires the resulting markovian dissipative 
time-evolution to be completely positive. 
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1. INTRODUCTION 



Recently, experiments based on neutron [1] and photon [2] interferometry have been 
proposed to test the hypothesis of noncontextuahty in quantum mechanics; the idea is 
to check whether a Bell-hke inequahty of the Clauser-Horne-Shimony-Holt form [3] is 
violated or not. Such an inequality is derived from the assumption that the measured 
values of physical observables are completely specified by the state of the system prior to 
measurement and that the actual measurement outcomes do not depend on the context, 
namely on whether other commuting observables are simultaneously measured. 

Differently from Bell-locality tests based on entangled physical systems, the above 
experiments involve two degrees of freedom of a same physical system; one degree of 
freedom is translational, related to the two possible paths followed by neutrons or photons 
inside the interferometer, the other is the spin (helicity) of neutrons (photons). 

Standard Bell-locality tests are not concerned with the time-evolution of the particles 
involved; only when there is lack of unitarity and loss of probability as in experiments 
based on decaying neutral K-mesons [4-6], the time-evolution becomes important. 

More in general, the dynamics is to be taken into account when the test particles 
behave as open systems S propagating through an environment £ to which they are cou- 
pled; in such cases, one usually traces away the degrees of freedom of £ ending up, under 
certain assumptions, with a one-parameter semigroup of linear maps Tf on the states of S 
represented by density matrices p. 

The maps Ft constitute a so-called reduced dynamics for the open quantum system 
S and embody the dissipative and mixing effects due to the environment S; they are not 
unitary, satisfy the forward in time composition law Tf^g = o for s,t > 0, and 
transform pure states into statistical mixtures. When they enjoy the property known as 
complete positivity they form a so-called "quantum dynamical semigroup" [7-10]. 

Physical consistency requires that the positivity of states of S, that is the positivity 
of the eigenvalues of the corresponding density matrices, be preserved for all times by 
any meaningful reduced dynamics F^; indeed, the statistical interpretation of quantum 
mechanics identifies these eigenvalues with probabilities. 

The property of complete positivity guarantees not only that the maps F^ preserve 
the positivity of the states of S, but also that the maps Ijv ® preserve the positivity of 
all states of the composite system Sn + S, for any A?"- level system, with In the identity 
operation on Sn [11,12]. Complete positivity of F^ is stronger than positivity and is 
intimately connected with quantum entanglement; indeed, positivity alone is not sufficient 
to ensure that Itv ® Fj preserve the positivity of entangled states of Sn + S. 

Noticeably, the standard quantum mechanical time-evolution generated by Hamilto- 
nian operators is unitary, reversible and completely positive. On the contrary, the physical 
literature abounds with dissipative, irreversible, reduced dynamics of quantum open sys- 
tems that are neither positive, nor completely positive, see e.g. [13-15]. 

In particular, in view of the abstract, experimentally uncontrollable coupling of the 
system of interest S with any A^-level system Sn, the argument that F^ should necessarily 
be completely positive may look as a mathematical convenience and a technical artifact, 
rather than a physical necessity [16]. In fact, the elimination of the environment degrees 
of freedom yields an equation of motion with memory terms that have to be eliminated 
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via suitable Markov approximations. It depends on how they are performed whether the 
resulting semigroups are physically consistent or not [17]. 

Recently, the issue of complete positivity has been reconsidered in the context of 
neutral meson dynamics, where a typical experimental situation is that of an entangled, 
singlet-like, state p of two K or B neutral mesons propagating according to the factorized 
time-evolution [18-21]. If the dynamical maps F^ are assumed to be not of the stan- 

dard Weisskopf-Wigner form, but modified by a noisy background of gravitational origin, 
it is showed that Ft has to be completely positive. Were it not so, physical inconsistencies 
as the production of negative probabilities, would affect F^ ® Ft[p]; moreover, these incon- 
sistencies cannot be dismissed as experimentally invisible because they might give rise to 
detectable effects [20] . 

In this paper we consider physical cases where the environment £ is given by a classical, 
fluctuating external field [13,22] and is a single open quantum system with two degrees of 
freedom. In such a context, the coupling is not between S and an abstract 2- level system, 
but between two degrees of freedom of S itself. Then, the physical meaning of complete 
positivity comes to the fore when we study a time-evolution of the form Ijv <E) F^. 

Namely, we study what happens if the interferometric apparatuses proposed for non- 
contextuality tests are placed in weak gaussian stochastic magnetic fields, or stochastic 
optical media, coupled to the spin, or helicity of neutrons, respectively photons, that will 
then propagate as open quantum systems. 

In particular, in the case of neutrons, we will consider in detail three possible choices 
of fluctuating magnetic backgrounds; these will give rise to reduced dynamics I2 ® F^ 
with Ti affecting the neutron spin degree of freedom covering all possible cases, namely 
to Vt completely positive, positivity preserving, but not completely positive and, finally, 
to F^ not even positivity preserving. The three possibilities depend on the properties of 
the stochastic magnetic field; it thus appears that by reproducing stochastic magnetic 
backgrounds with the qualities of the three cases we are referring to, one would be able to 
experimentally study the characteristics of the various reduced dynamics. 

Further, we study how the Clauser-Horne-Shimony-Holt inequality is modified by the 
stochastic magnetic field and show that complete positivity of the time-evolution Ijv ®Ft is 
necessary for a consistent physical description; otherwise unacceptable negative eigenvalues 
appear in the time-evolving physical states describing entangled degrees of freedom. 



2. ENTANGLEMENT AND NON-CONTEXTUALITY TESTS 

In the following we shall refer to noncontextuality tests using neutron interferometry 
[1] , photons involving similar arguments [2] . In the experimental setup proposed in [1] , an 
incoming beam of neutrons with spin along the positive ^-direction passes through a beam 
splitter with transmission and refiection coefficients p and q with |pp + |qP = 1. The beam 
gets divided into two components that follow two spatially separated paths u and d. 

Both the spin and the translational degree of freedom are described by 2-dimensional 
Hilbert spaces, the former with basis vectors ] t«) and | J.^), the latter with basis vectors 
\il)u) and iV'd) corresponding to the two possible macroscopic paths. Making the beam 
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tt-component undergo a spin-flip | f^) i— > | J,^), an initial beam state is prepared, 

m^p\i^u)®\iz) + q\i^d)®\U) , (2.1) 

which then propagates inside the interferometer. A vector state as above corresponds to 
the one-dimensional projector := |^')(^'|, 

= \pfPi (g) Q2 + \q\'^P2 ® Qi + pq* Ps^Qa + P*qPA ® Qs , (2.2) 

with 

Pi ■■= \i^u){i^u\ , P2 ■■= MM , P3 ■■= \i^u)M , P4 ■■= Mii^ul , (2.3a) 

Qi := I I , Q2 := I iz){iz I , Q3 := I h){iz \ , Qa := \ iz){h I • (2.3b) 

More in general, neutron beam states are not pure as p.^, rather statistical mixtures de- 
scribed by density matrices 

4 

P(2) — Yl ® ' (2.4) 

that is by 4 X 4 hermitian, normalized, positive matrices, whose positive eigenvalues sum 
up to one (Trp(2) = 1). Positivity is crucial for the statistical interpretation of quantum 
mechanics where the eigenvalues of density matrices play the role of probabilities. 

Neutrons spend a typical time t within the interferometer during which they may be 
subjected to external influences resulting in a dynamical change p(^2) ' — * P{2){t) of their 
state. At the exit of the interferometer, a second beam splitter recombines the translational 
components and shifts the u component by an angle 



IV'd) J ^ \\M^,v)) J le-^'^cosi? -smi&J Ui/jd) 



(2.5) 



with reflection and transmission probabilities cos^ -d and sin^ -d. Consequently, the neutron 
beam state emerging from the interferometer is 

p ^ {Uii}, ^) ® I2) P(2)(t) {U*{i}, ^) ® I2) , (2.6) 

where U{i}, </?) is the unitary matrix in (2.5), U*{i}, 99) its adjoint and I2 is the 2x2 identity 
matrix. The two components of the exiting beam are then intercepted by two counters Cu,d 
plus spin-analyzers ^n,-!! that record how many neutrons reach them with spins polarized 
along suitable directions ±n = ±(ni,n2,n3) in space. The frequencies of counts give the 
expectations 

Oi'^^i^^if) :=Tv(p(2)(t)P,(7?,(^)®Qn) , (2.7) 

where Pj{'&, (p) U*{'i}, ip) Pj U{'i}, ip), Pj, j — 1, 2, are as in (2.3a) and represent observ- 
ables that are chosen by setting the angles (p characteristic of the beam splitter. Further, 
Qn '■= I Tn)(Tn | projccts outo a statc with spin along the direction n. 
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Since translational and spin observables commute, the observables 



A{^,ip):=Pi{^,ip) - P2{^,^), S(ii) := - Q_„ (2.8) 

also commute and have eigenvalues ±1. Choosing angles '!?i,2, <fi,2 and polarization direc- 
tions ni^2, one constructs commuting observables Ai, Bj, i,j = 1,2, called dichotomic [1]. 

In the hypothesis of noncontextuality, the possible outcomes ±1 of a measurement 
of Ai, respectively Bj, are predetermined by the state P{2)it) independently on whether 
Sj, respectively Ai, are simultaneously measured with Ai, respectively Bj. From such 
assumptions a Clauser-Horne-Shimony-Holt inequality can be derived for the mean values 



Ct{d, ^; n) := l¥(p(2) ^) ® B{n)) 



(2.9) 



with four possible configurations of the control parameters i?, ip and n, 

Ct(^?i, Vi\ ni) + Ct(^?i, <^i; n2) + Ct{^2. ni) - Ct(i?2, V2\ ns) | < 2 . (2.10) 

From (2.9), it turns out that the quantities Ct{'&, ip; n) can be measured by frequencies at 
counters plus spin-analyzers iCu.Sn) and (Cd,<S'_n) with the beam splitter at the exit of 
the interferometer set at angles d and ip. 

If inequality (2.10) is violated, the hypothesis of noncontextuality upon which it was 
derived, cannot hold. Whether it is so or not can be checked in highly efficient experiments 
where the entanglement between translational and magnetic degrees of freedom is exploited 
[1] (for a similar argument involving photons see [2]). 

Interestingly, by setting appropriately the angles d and (p of the beam-splitter and 
the polarization direction n of the spin- analyzers, also the entries of the state P(2)(^) = 
Y^i j Pij{t)Pi ® Qj, can be measured. From the entries, one has access to the eigenvalues 
of the beam state after travelling through the interferometer and thus to the effects of the 
time-evolution inside it. Indeed, from (2.3, 2.4) and (2.6) it readily follows that 

Piiit) = 0^^(0,0) , pi2(t) = 0^-^(0,0) , p2i(t) = Ot''(0,0) , p22(t) = c?''"'(o,o) . 

(2.11) 

Furthermore, since operators Pi®Qj with i,j = 3,4 are not self-adjoint, their expectations 
can only be measured indirectly, through the mean values of the projectors 

P± := PiA^. 0) , P±. := Pi,2(^, -\) , (2.12a) 
^ . I Tz) ± I i^) (T^ I ± I r, . I T^) I i^) (T^ I I 

Q±..= ^^l 71^'^±-= — T2 vi — • ^ ^ 



Then, from 

~ , = "^-^ '-^-^ = Ql , (2.13) 



P+ - P- +iP+i-iP-i ^ Qx-Q-x+iQy-iQ-y 



5 



one obtains expressions for all other entries p-j^ (t) := Tr(^p(2)(t)P* ®Q'jkj- We quote one 
of them which will be needed in the sequel ( the others are reported in Appendix A), 

P43(t) = Z + I Z 



0''"(f , 0) - o^-^(f , 0) . o,^'^(f , 0) - o'r'il, 0) 

2 2 ^ 2 2 



+ ^ 



, '-^t ^4' 2/ '^t l4' 2/ '^t ^4' 2/ '-^t '^l' 2/ 



4 4 
^^2,x/n _7r\ _ /n2,— aj/TT _7r\ (r\^,y ( -k _ tt ^ _ /-/-)2,— j/ / tt _ tt \ 
'-^t ^4' 2/ '-^t V4' 2/ I '-^t ^^4' 2/ '-^t V4' 2' 



-i ' . ^ ^ + 



(2.14) 



3. OPEN SYSTEM DYNAMICS INSIDE THE INTERFEROMETER 

Neutron interferometry has proved an extremely powerful tool to investigate gravi- 
tational, inertial and phase-shifting effects occurring inside the interferometer [23-29]. In 
the following we will show that neutron interferometry might also be used to investigate 
the notion of completely positive open system dynamics. In order to do that, we consider 
the case in which neutrons while propagating inside the interferometric apparatus are sub- 
jected to weak time-dependent, stochastic magnetic fields coupled to their spin degree of 
freedom. 

We assume the time-dependent Liouville-Von Neumann evolution equation for the 
2x2 density matrix S describing the spin degree of freedom to be of the form 

atE(t)-(Lo + Lt)[E(t)] , 

Lo[S(t)]:=-i [5^(73, , L4S(t)]:=-i[v(t)-5=,S(t)] , ^^"^^ 

where a = (cri, cr2, 0-3) is the vector of Pauli matrices, V(t) = (V\(t), V2(t), V3(t)) is pro- 

portional to gaussian stochastic magnetic field and Hq := — 0-3 is due to the coupling to 

a static magnetic field along the 2;-direction. Furthermore, we assume V(t) to have zero 
mean, (V(t)) = 0, and stationary, real, positive-definite covariance matrix W(t) = [Wjj(t)] 
with entries 

W,,{t -s) = mmi^)) = W*,{t -s) = W,,{s - t) . (3.2) 

Because of the stochastic field V{t), the solution E(t) of (3.1) is also stochastic; an effective 
spin density matrix p{t) :— (S(t)) is obtained by averaging over the noise. At time t = we 
may suppose spin and noise to decouple so that the initial state is p :— (S(0)) = S(0). In 
order to derive an effective time-evolution for p{t), we follow the so-called convolutionless 
approach developed in [13]. 

We average over the noise in the interaction representation, where we set 

S(t) := exp (-tLo)[E(t)] = e"**^" E(t) e**^° , (3.3a) 
p(t):=(E(t)) and Lt := e'*^" Lte*^° . (3.3b) 
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The result is 



~p{t) = Y,M2k{t)[p\ 



where 



fe=o 



Mk{t)[p\:^ I dsi [ ds2--- [ dsk {L{si)L{s2) ■ ■ ■ L{sk))[p] 











(3.4) 



(3.5) 







Only even terms contribute to (3.4) because the stochastic field is assumed to be gaussian. 
Denoting by Mt the formal sum in (3.4), a resummation gives 



= ( M2it) + M^jt) - M2{t)M2{t\ 



2nd order 



4th order 



+ ^ )W)] 

higher orders 



(3.6) 



Since the action of the magnetic field on the travelling neutrons is, by hypothesis, weak, one 
can focus on the dominant first term in the expansion, neglecting higher order contributions 
[8-10]. By means of (3.2), the second order contribution can be worked out explicitly. 



Mt)W)] = - E f ^sWabis) k(t), k(t - s) , p(t)ll . (3.7) 

Returning to the Schrodinger representation and using the statistical independence of the 
Hamiltonian Hq from the stochastic field, it follows that p{t) := exp(tLo)[p(t)] solves 



where 



dtp{t) ^ -i Ho , p{t) - E Cijit) Oi, Oj , pit) 



3 .t 



cosojot — sintuoi 
U{t) I sinc^o^ cosLOot 
1 



(3.8a) 
(3.8b) 



(3.9) 



is the unitary matrix U{t) = [Uij{t)] such that e *^°[a-j] = ^ Uij{t)aj 



i=i 



From (3.2) and (3.9) it follows that C{t) := [Cij{t)] is a real matrix and can thus be 
decomposed into symmetric and antisymmetric components. Correspondingly, the second 
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term on the right hand side of (3.8a) sphts into a commutator with a Hamiltonian and a 
purely dissipative contribution 



J2 Cij{t)[a,, [aj , p(t)]] =i[Yl C$it)eijk(Tk , pit)] (3.10a) 

3 ^ 

E 2Cg(t)(-{a.a, , p{t)] - <J,p{t)<j^ , (3.10b) 



where C;J'^(t) = (C,,(t)TC,,(t))/2. 



When the couphng between system and stochastic field is weak, the memory effects 
in (3.8a) should not be physically relevant; therefore, the use of a Markov approximation 
is in general justified.^ In practice this is done by extending to +oo the upper limit of the 
integral in (3.8b); the resulting equation of motion has no explicit time-dependence, 

dtp{t) = -i[H , p(t)] + LdW)] , (3.11a) 

3 

H = Hq + Hd, Hd E ^ij^ijkfyk (3.11b) 

i,i,fe=l 

3 ^ 

LdW)] := E ' + ^jP{t)^) ' (3-llc) 

where C,4 and L^, are the entries of the real matrices 

:= r , (3.12a) 

(3.12b) 







Cd := ds {W{s) U{-s) + U{s) W(-s)) 



The Hamiltonian contribution is skew-symmetric, while the purely dissipative one, 
Ld[-], is symmetric; this latter makes the time-evolution irreversible, but preserves prob- 
ability because Tr^LD[p(t)]^ = 0. The solutions of (3.11a) thus constitute a semigroup of 

linear maps Vt : p i — > p(t) :— Tt[p\ such that r^+g = o F^, s, t > and Tr(Ft[p]) = Trp. 
It remains to be checked whether they preserve the positivity of spin density matrices, 
that is whether the F^'s are, in short, positive maps on the spin states, which is the first 
request for physical consistency. 



t More precisely, one can show that a linear, local in time subdynamics is the result 
of a limiting procedure in which the coupling constant g between system and external 
environment, and the ratio r/T between the typical time scale of the system and the 
decay time of the correlations in the environment, become small. [8-10] The quantities g 
and r/T regulate both the expansion (3.6) and the Markovian approximation of (3.8b). 
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There is, however, a further constraint that has to be respected for physical consis- 
tency. Indeed, if the neutron interferometer is placed in a stochastic classical magnetic 
field of the kind described above, the translational degree of freedom is not affected and 
the effective state P{2){t) at the exit from the interferometer will be P{2){t) — (I2 ®rt)[p(2)], 
where I2 denotes the identity operation on the first factors in (2.4). It turns out that the 
positivity of the maps Tt does not guarantee the positivity and thus the physical consis- 
tency of the maps I2 ® F^; for this the stronger notion of complete positivity has to be 
imposed on the maps Ft- 

We shall later investigate these notions in more technical detail, for the moment we 
observe that generating fluctuating magnetic backgrounds with certain decaying proper- 
ties of their covariance matrix (3.2), one may have experimental access to some physical 
situations of theoretical interest that we will present below. 



3.1 White noise 

The stochastic magnetic field has white-noise correlations 



Wij{t -s) = {Vi{t)Vj{s)) = Wij S{t - s) , (3.13) 

where W := [Wij] is time-independent, symmetric and positive-definite. Then, Ca = and 
jCd = W; furthermore, writing W — with A = [ay] real and symmetric, the dissipative 
term in (3.11a) reads 

L^[p]=5^AfcpAfe - l{Y,Al,p] , (3.14) 

k k 

with self-adjoint Ak := Xli=i '^kiO'i- This is a particular instance of Lindblad's theorem 
[30,31] which states that a family of linear transformations Ft : p ^ p{t) on the D- 
dimensional density matrices is a quantum dynamical semigroup of probability-preserving, 
completely positive maps if and only if it is generated by the equation of motion 

dtp{t) = -i[H, p{t)] - ^[j2AUe,p{t)] + J2Aep{t)A} , (3.15) 

£ £ 

where the Ag^s are D-dimensional matrices with adjoint A| such that the series are (norm-) 
convergent. On the other hand, if Ld in (3.11a) is as in (3.15), then the corresponding 
jCd is positive definite. Further, note that, because of (3.13), all higher order terms in 
the expansion (3.6) identically vanish, so that the evolution equation (3.15) is in this case 
exact [22]. 

3.2 Diagonal covariance matrix 

The stochastic field has no off-diagonal correlations, while 

mt)v^is)) = mt)V2is)) = ^^SiV^i*-^! , mmis)) = ^^sie-^i^-^i , (3.16) 
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where are constant magnetic field intensities and g is proportional to the neutron 
magnetic moment. Then, 



Ca 



A2 + a;2 




-D 



Setting 



Aa; := 



A2 + 



a;? 



7 := 



V 



a := 





XB'l 





\ 





(3.17) 



A2 + 



a;, 



(3.18) 







the matrix in (3.12b) becomes 




The reason for such a parametrization will become clear in the next section 

>l(t) P3(t) 



(3.19) 



Given the entries of the spin matrix p(t) 
satisfy the Bloch-Redfield equations [32] 



are readily showed to 



Pi = -7Pi + 7P2 , P3 = -*(<^o + Aa;)p3 - 2ap3 



(3.20) 



and p2 = —pit Pa = (ps)* ■ The coefficients 27 and 2a are the inverse of the relaxation times 
Ti and T2 of the diagonal, respectively off-diagonal elements of p{t); from the positivity 
condition 2a — 7 > 0, it follows that: I/T2 > l/2Ti. In [13,14] it is showed that this typical 
order relation can be reversed by setting ^3 = and keeping 4-th order terms in (3.6). 
In such a case, however, I/T2 < l/2Ti implies a — 7/2 < and jCd in (3.19) is no longer 
positive-definite. 

By Lindblad's theorem, the argument of case 1. implies that the corresponding dynamical 
maps generated through (3.11a) can not be completely positive. We shall see in the 
next section that at least they preserve positivity. 



3.3 Single component field correlation 

The stochastic magnetic field is along the x-direction, V{t) = (Vi(t),0,0), with 



mt)Vi{s))^g'B'e 



2 d2 „-A|t-s| 



(3.21) 



Then, 



Ca 



9^u;oB^ 




A2 + 



2A ujQ 
c^o 




(3.22) 
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Unless ojq = the matrix >Cd is not positive-definite. By analogy with the parametrization 
of the previous example, we set 

Then, the corresponding Bloch-Redfield equations for the entries of p{t) read 

Pi = -IPi + 1P2 , h = -i{(^o + Aa;)p3 - 7P3 + 7P4 + 2ibp4 , (3.24) 
and p2 = -pi, P4 = (ps)*- 



4. COMPLETE POSITIVITY VS SIMPLE POSITIVITY 

As already remarked, physical consistency demands that Ft preserve the positivity of 
initial density matrices p describing the neutron spin degree of freedom. 

In order to check whether this is so in the preceeding cases, it is convenient to decom- 
pose the spin density matrices p by means of the Pauli matrices aj, j = 1,2,3 plus the 
2-dimensional identity matrix ctq, namely p = S^=o P'^'^m- such a way, density matrices 
can be represented as 4-dimensional ket-vectors |p) = {p^ , , , p^) , where 

Pa - P3 „3 Pi - P2 



P° = ^^, P^ 



P3 + P4 2 

rt ' P 



2i 



(4.1) 



Since they operate linearly, the commutator and the purely dissipative term Ld[ 
in (3.11a) act on the vectors |p) as a skew-symmetric matrix 



n = -2 



























y 



(4.2) 



respectively as a real, symmetric matrix 



V = -2 



^0 ON^ 

a b c 

b a P 

VO c 7/ 



(4.3) 



The connection with in (3.12b) is readily derived to be 



a + 7 — a —2b —2c 
= - I —2b a + 'f — a —2j3 

-2c -2(3 a -Ha -7 



(4.4) 



which explains the parametrization used in the previous section. 
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In this representation, the time-evolution equation (3.11a) reads 



dMt))^{n + vMt)) . (4.5) 

In order to find necessary and sufficient conditions for Tf to be positivity preserving, we 
now proceed in a few steps. 

First, since = 0, the trace is conserved, thus p^{t) = 1/2; therefore, the positivity 
of Tt[p] is ensured if Bet[p{t)] = 1/4 - E5=i(p-^)^ > 0. 

Second, the time-derivative of the determinant at i = must be positive whenever 
Yl^=iiP^)'^ = 1/4 (so that Det[p] = 0), otherwise one eigenvalue would become negative 
for t > 0. Using (4.3,4.5), we thus get the necessary condition 

= _2 ± P,,,v- > . (4.6) 

By varying p^ , while keeping J^jiP'')'^ — 1/4, it follows that 

a b c\ 

V^^^ := -2 I 6 a /3 (4.7) 

c (3 7/ 

must be negative definite which in turn implies 

a > , aa > b^ , Detp(^) > . (4.8) 

Third, conditions (4.8) arc also sufficient for Tf to preserve positivity. In fact, since 
—T> > 0, we can write —V = with B symmetric. Then, the term in the right hand side 
of the equality in (4.6) is given by ||i3|p)|p. Let us suppose Det[p(t')] < 0, at time t' > 0; 
it follows that Det[p(t*)] = at some time t* such that < t* < t'. Thus, B\p{t*)) = 0, 
otherwise Bet[p{t)] > for t >t*; but this implies \p{t)) = \p{t*)) for all t>t* under the 
time-evolution 

\p)^\p{t)) = e'^\p) = J2t^''\p) ■ (4.9) 

Therefore, as well as the standard dynamics generated by a Hamiltonian operator, the 
dynamics (4.9) is positivity-preserving; via the Lie- Trotter product formula, it then follows 
that the time-evolution generated by Ti -\-V^ 

Gt := exp (t{n + P)) = Jirn (e^/^^e*/^^)"" , (4.10) 

preserves positivity, too. 

However, even if Tt preserves positivity of the states describing the neutron spin degree 
of freedom, this does not guarantee that I2 ® T^ preserves the positivity of states in which 
the spin is entangled with another degree of freedom. For neutrons in the interferometric 
apparatus of the previous section, such a request is crucial since the maps I2 ® Tt tell us 
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how do evolve in time states p(2) describing both the magnetic and translational degrees 
of freedom. A theorem of Choi [12] states I2 <S> to be positivity preserving if and only if 
is completely positive. 
Among the neutron states p(2) propagating through the interferometer, those without 
correlations between translational and spin degrees of freedom are of the form Pspace ® Pspin 
or are linear mixtures of them. If the Tt are positivity preserving these tensor-product 
states remain positive in the course of time; indeed, 

< Pspace ® Pspin ' > h ® T^t[Pspa,ce <8 Pspin] = Pspace <8 T^PspiJ > . (4.11) 

However, it is not so for entangled states. 

Let us take p — —q — l/\/2 in (2.2), so that the initial beam state is antisymmetric 
in the two degrees of freedom. If we ask I2 (8) Tt to preserve the positivity of p^, it must 
hold that 

A^{t) := ($|(l2 (^r^) [p^]|$) > (4.12) 
for all If $ is orthogonal to the fact that A$(0) = implies 

= ($|(l2 0Li^)[p*]|$) >0 . (4.13) 

t=o 

By varying $ in the 3-dimensional subspace orthogonal to ^' we obtain the inequalities 



2R = 


q: + 7 - 


a > , 


RS > 


2S = 


a + 7 - 


a > , 


RT > 


2T = 


a + a — 


7>0, 


ST>p^ 



(4.14) 



RST >2bcP + R/3^ + Sc^ + Th'^ . 

These inequalities are stronger than the ones in (4.8) and must necessarily be satisfied 
if we want to avoid that the maps I2 <8) Tt become physically inconsistent by generating 
negative probabilities out of initially entangled beam states. 

Furthermore, using (4.4), inequalities (4.14) amount to the positivity of the matrix 
Cd = [Lij] of the coefficients of the dissipative term in (3.11c) and thus they imply the 
complete positivity of the time-evolution Vt- 

Let us now discuss the three cases introduced in Section 3. 

4.1 White noise 

As already noticed, the evolution equation (3.11) is now exact. The matrix T)^^^ takes 
the general form (4.7), so that provided the inequalities (4.14) are satisfied, the integrated 
time evolution Tt results completely positive. 

4.2 Diagonal covariance matrix 

In this case one finds: 

ja 0\ 

= -2 a . (4.15) 
Vo 7/ 
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The corresponding dynamics is completely positive only when I/T2 > l/2Ti (a > 7/2); it 
is positive, but not completely positive, when I/T2 < l/2Ti (a < 7/2). 

Analytic solutions of the equation of motion corresponding to (4.15) are readily cal- 
culated; in vectorial representation one has 

p'{t) = P° 

pi(t) = e~^°-*\ cosujt - p'^sinujtX 

(4.16) 

p'^(t) = e j sin Lot + p^ cos Lot \ 
p\t) = e-'-r^p' 



with Lo :— loq + Au). The difference between complete positivity and positivity shows up 
in different order relations between the decay diagonal and off-diagonal relaxation, that is 
either a > 7/2 or a < 7/2. 

However, the true physical meaning of complete positivity of Tt becomes evident 
when the state in (2.2), with p — —q = l/\/2, evolves in time according to I2 ® F^. In 
vectorial representation \Qi) = 1/2(1,0,0,1), IQ2) = 1/2(1,0,0,-1), IQ3) = 1/2(0, l,i,0) 
and \Q4) = 1/2(0, 1, — i,0), thus, using (4.16) one obtains 

P(2){t) = {h<»rt) [p^] = ^(Pi(8)Q2(t) + P2®Qi{t) - Pz®QA{t) - PA®Qz{t)) , (4.17) 



with 



Qi{t) 







1 - e-27* 



Q2{t) = 7. 



1 / 1 _ e-27* 



1 + e- 



-2-yt 



1 






1 



(4.18) 



Therefore, 



P{2)it) 



fE_{t) \ 

E^{t) F{t) 

F*{t) E+{t) 

V E_{t)J 



E± := 



1 ± e-27* 



F := 



Q—t{2a—iui) 



(4.19) 

The eigenvalues of the state P{2){t) at the exit of the interferometer are Ai^2(^) = E_{t) 
and 

A±(t) = . (4.20) 

Let a < 7/2, that is let to be positive, but not completely positive; then, since A_(0) = 
and dA+(0)/dt = (2a — 7)/2, there is a whole range of t where X+{t) < and P{2){t) loses 
physical meaning. On the other hand, if a > 7/2, X+{t) > 0, for all t. 
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4.3 Single component field correlation 

The matrix V^^^ in (4.7) has now also ofF-diagonal terms: 

P(=^) = -2U 7 0, 6=-^; (4.21) 
VO 7/ 

only when these are zero, i.e. Aoj = 0, positivity is preserved. This can be seen also by 
considering the integrated time evolution; in the vectorial representation, one explicitly 
finds: 

p\t) =e-^^\p^( cosh St + ^sinh^t) - '^i±^smhSt} 

_nr. (4-22) 

p^(t) =e-^*|p^ ^^:^^sinh5t + p^ (cosh St - |sinh5t)} 

with S := y^7^ + (46^ — w^). One notes that, besides developing negative eigenvalues 
because of lack of positivity, this evolution cause p^^''^\t) to diverge with t i-^ +oo, when 
S > J. 

Although apparently formal, these results are far from being academic: indeed, as 
already stressed at the end of Section 2, the entries of p(2) are directly accessible to the 
experiment.^ By modulating a background magnetic field close to the stochastic properties 
investigated in the previous three cases, one might reproduce experimentally the conditions 
for three different reduced dynamics and check their consequences. 

Then, one may conclude that reduced, markovian time-evolutions Tt must be not 
only positive, but also completely positive, since lack of any of these constraints results in 
experimentally detectable inconsistencies. 

Clearly, the use of one reduced dynamics instead of another depends on the markovian 
approximation used to derive it and whether, given the properties of the stochastic field, 
it was justified or not. It thus seems appropriate to conclude that, whenever a semigroup 
composition law is expected, the physically appropriate markovian approximations are 
those which lead to reduced dynamics consisting of completely positive maps [17]. 



5. NONCONTEXTUALITY AND DISSIPATION 

We now examine to what extent the Clauser-Horne-Shimony-Holt inequalities (2.10) 
are modified by the presence of a stochastic magnetic field with covariance matrix as in 
the cases studied in Section 3. We first notice that the mean values (2.9) can be written 

C,(79,^:n) = 2(o^"(^,(^)-0^-"(^,(^)) - IV^ (pspin(t)S(n)) , (5.1) 

T For instance, in the case 4.2 above, it follows from (2.11) that E_{t) = Ol'''{0,0) and 
E+{t) = of''' {0,0), while F{t) coincides with the expression in (2.14). 
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where Tr2 denotes the trace over the spin degree of freedom and Pspin(^) = T^iP{2)it), T^i 
denoting the trace over the translational degree of freedom. 

For sake of simphcity, we again consider an initial beam state with p = —q = l/\/2; 
then, 

^; n) = 2 ^) - v)) 

= sin^^ lV2(g2(t)S(n)) + cos'^^* IV2 (gi(t)S(n)) 

Further, wc shall take — — and h'^ — u)q/2 in the hamiltonian part (4.2) of the 
evolution equation (4.5). 

From (4.10) it follows that, in vectorial representation, the time-evolution operator 
acts on initial states \p) as 



^1 ^^ 

Quit) g,2{t) gi3{t) 

g2l{t) Q22{t) Q2z{t) 

Vo g^iit) gs2it) gzz{t)I 



(5.3) 



Also, the mean value of an observable X = X];/=o -^^^i' with respect to p is given by 
Tr(Xp) = ^Y^^^qX^ p'^ ■ Further, the observable B{n) — — Q-n relative to n = 
(ni,n2,n3) corresponds to the vector \B{n)) = (0, ni, n2, na); thus one computes 



Tr 



2 (Qi,2(t)5(n)) = ±G{t) • n , Tr^ {Q^{t)B{n)^ = F{t) 



n 



where 



G(t) = (gi3{t), fe(t), gssit)) and 

F{t) = {gii{t)-igi2{t), g2i{t)-ig22{t), gsiit) - ig32{t)) 

Finally, the mean values (5.2) read 

Ct{'d,(p;n) = n- cos2i?G(t) - sm2^ne(e-''^F{t) 



(5.4) 

(5.5a) 
(5.5b) 



(5.6) 



Wc shall now discuss the explicit behaviour of Ct{'d, (p; n) in the three cases introduced in 
Section 3, and further analyzed in the previous section. 

5.1 White noise 

Though analytic expressions of are obtainable in the general case of a stochastic mag- 
netic field with white noise correlations, these are rather involved and scarcely illuminating. 
More conveniently, one may suppose the dissipation T> in (4.3) to be small in comparison 
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to the Hamiltonian contribution due to Hq. In practice, we assume the parameters a, b, 
c, a, (3 and 7 small with respect to w = cuq + ~ and proceed with a perturbative ex- 
pansion (for more details compare [18,19,21]). To first order in the parameters the entries 
of Qt are as reported in the Appendix B. Using them, one calculates 



^ m 4|C| . ivot fujot 
Gi(t) = -^sm — cos (— 

4|C| . ujot . /u>ot 

— — - sin sm I 

G3(t) = e-2^* , 



G2(t) = ^sinr^sin(^::^ + 0c) , 



where |Cp = + and tan 0c = P/c, 



/ \ \B\ 

7^e(e-*'^Fl(^)) = e-^''+'')* cos(a;ot - ^) + ^sina;ot cos((^ + 0b) 



B\ 



ne(e-''^F2{t)) = e-("+")* sin(a;ot - ^) - — sinc^o^ sin(^ - 0b) 



7^efe-^^F3(^)) =-^sin^tcosf^t-^ 



where = (a — ct)^ + 46^ and tan 0s = 2b/ {a — a). It thus follows that 

Cti'd, (p; n) = ni — cos 2'?? sin ^tcos (-^t + 0^) 
L o^n 2 V 2 / 



sin2T? (e"*'^''"'""^ cos(a;o^ - V') + — - sina;o* cos(<^ + 0b) ) 

■4|C| ^ a ■ '^0 . /'^O , \ 

COS zv sm —tsm[—t + 0c) 

. ujQ 2 2 

-sin2i?(e~*'^''+'') sin(u;ot - (fi) - - — - sin wot sin((/p - 0b) ) 
+ ns e-2^*cos2i$' "^1^' 



+ n2 



sm 2v sm — t cos (^"2" ~ f ~ 9c j ■ 



(5.7a) 

(5.7b) 
(5.7c) 



(5.8a) 
(5.8b) 
(5.8c) 



(5.9) 



5.2 Diagonal covariance matrix 

For the stochastic magnetic fields with correlation matrices as in (3.16) we have 

G{t) = (0, 0, e-2^*) , F{t) = e-*(2«-^'^> (1, -i, 0) . (5.10) 
With these expressions one easily derives 

Cti-d, (fi; n) = e"^^* cos 2t9 nz - e'^"* sin 2-3 (ni cos{ujt - (fi) + n2 sm{ujt - (fi)^ . (5.11) 
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5.3 Single component field correlation 

The computation is similar in the case of magnetic field with covariance as in (3.21). One 
finds: 

G{t) = (o, 0, e-^T*) , (5.12a) 
F(t) = e-'T* (^{coshSt + ijsmhSt) (1, -i, 0) + ^^^^sinh^t (1, i, 0)) ,(5.12b) 

from which one esaily obtains, assuming S > 0: 
Ct(^,ip;n) = e-^^* cos 2^713 



e"^* sin 2i? 



rii ^(cosh 6t + sinh 6t) cos (p + '^^^-^ — - sinh St sin ip^ ^3) 
+ ^^-2 ^( — cosh + sinh St) sin (p + ^ sinh St cos 

The lack of positivity preservation which characterizes the time-evolution leading to (5.13) 
manifests itself in that the quantities Ct{p., n) diverge with large t when 5 > 7. 

Expressions (5.9), (5.11) and (5.13) agree with those used in [1,2] when there is no 
dissipation, namely putting a = 6 = c= Q! = /? = 7 = 0, 

(^^(■i?, (/?; n) = — 77-1 sin2'j?cos(a;ot — v^) — 77-2 sin 2'i?sin(a;ot — (y?) + n3COs2'i? . (5.14) 

Notice that the unitary time-evolution generated by the Hamiltonian Hq contribute to a 
time varying redefinition of the angle 

Concerning the issue of complete positivity vs simple positivity, in expressions (5.9) 
and (5.11) the two possibilities manifest themselves in different relaxation properties due to 
whether inequalities (4.14) or (4.8) are fulfilled. No physical inconsistencies may affect the 
mean values Ct(^?, n); indeed, negative probabilities may result in negative mean values 
of positive observables only if the latter are entangled. In the case of the quantities involved 
in inequality (2.10), the observables are factorized, P\^2{^-,'p) ® Qn and the positivity of 
their mean values is preserved even when Ft is only positive and not completely positive. 

This can be seen as follows. To the Schrodinger time-evolution P{2)it) — i}-2®^t) [P(2)]; 
there corresponds the Heisenberg time-evolution of observables X(2)(t) = (I2 ® ^t) [-^(2)]) 



Tr 



((I2 ® r^) [p(2)] ^(2)) = 1V(P(2)(I2 ® ^t) [^(2)]) • (5.15) 



The maps O^, dual to F^, form a semigroup of dynamical maps that transform positive 
observables into positive observables, if the Ft's preserve the positivity of states. Conse- 
quently, even when the initial state p{2) is entangled and the Ft's positivity preserving, but 
not completely positive, it turns out that 

Tl(p^2){t)Pl,2{'&,V)®Qr) =Tl(p^2)Pj{'^,V)®^t[QrS) >0 . (5.16) 
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6. CONCLUSIONS 



Complete positivity is a property of quantum time-evolutions which is enjoyed by 
the standard dynamics of closed quantum systems generated by Hamiltonian operators, 
but not automatically by the more general reduced dynamics describing time-evolution 
of open quantum systems in interaction with suitable environments. Complete positivity 
is intimately related to the phenomenon of quantum entanglement, between two different 
systems, but also between two different degrees of freedom of a same physical system. 

In this paper we have considered the two entangled degrees of freedom, translational 
and rotational, of a beam of neutrons travelling through an interferometric apparatus de- 
voted to checking the hypothesis of noncontextuality. We have studied the consequences 
of placing the interferometer in a stochastic, gaussian magnetic field weakly coupled to 
the spin degree of freedom that provides an experimentally controllable environment. As 
explained in Section 3, the same Markov approximation naively yields a semigroup of dy- 
namical maps I2 , where only the spin degree of freedom evolves in time; by varying the 
decay properties of the external field correlations, these maps turn out to be alternatively 
completely positive, simply positivity preserving, not even positivity preserving. 

The noncontextuality tests proposed in [1,2] are based on the Clauser-Horne-Shimony- 
Holt inequality (2.10) without time-dependence, that is with t = 0. The presence of a 
fluctuating magnetic field induce relaxation on the spin degree of freedom with strength 
and properties depending on those of the field. Typically, the mean values in the inequality 
are damped and make it more difficult to be violated. However, in presence of stochastic 
fields yielding reduced dynamics that do not preserve positivity, the inequality might be 
dramatically violated because of possible mean values diverging in time. 

This latter possibility is a manifestation of the fact that any physically consistent 
time-evolution Ft must preserve the positivity of spin states in order that the eigenvalues 
of the corresponding spin density matrices might at any time be used as probabilities, in 
agreement with the statistical interpretation of quantum mechanics. If the F^'s preserve 
the trace of spin density matrices, but not their positivity, spin states may evolve in time 
in such a way that some of their eigenvalues become negative, while others greater than 
1, without upper bounds. It is this physically unacceptable phenomenon that leads to 
diverging mean values. 

The request of positivity preservation by the maps with respect to spin states is 
thus unexcapable, but it is not enough to avoid physical inconsistencies when the time- 
evolution maps I2 <S> Tt act on states p(2) with correlations between spin and translational 
degrees of freedom. 

Inequality (2.10) does reveal the difference between completely positive and simply 
positivity preserving F^, but only as long as the relaxation characteristic are concerned, 
without any further effect (as the divergence of some contributions to the inequality). In 
fact, the positive observables in (2.10) arc factorized, that is they incorporate no entan- 
glement between the translational and spin degree of freedom. Even if the initial state 
does incorporate entanglement, it nevertheless follows that the mean values of factorized 
observables remain positive and bounded. 

However, the interferometric apparatus proposed in [1,2], might also be used to mea- 
sure the entries of the states of the neutron beam at the exit of the interferometer. In 



19 



this way, one might have access to the spectrum of an initially entangled state after being 
subjected to the effects of the stochastic magnetic field. 

In the case of fluctuating flelds yielding reduced dynamics that preserve positivity, 
but are not completely positive, the theoretical predictions indicate the appearance of 
negative eigenvalues, that is of negative probabilities, in the spectrum of the entangled 
exiting state. The fact that they are, in line of principle, detectable experimentally, does 
not allow to dismiss such an occurrence as practically negligible. Rather, it forces to 
reconsider the Markov approximation used to derive the time-evolution and to select as 
physically consistent only those providing completely positive reduced dynamics. 
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APPENDIX A 



Pilit) 

P^3{t) 



0^-(O,Q) - 0^-^(0,0) ^ . Ol'-' (0,0) - Ol'\0,0) 

—I— Z 



(!)?'^(Q,Q) -0^'-^(O,Q) ^ .0,^'-^(O,O) -0?'^(O,Q) 

_ —I— z _ 



o;'^(f , 0) - o?'^(f , 0) o;'^(f , -f ) - o?'^(f , -f : 



0,^'--(f,Q)-0,-'--(f,Q) ^ 0,^'--(f,-|)-0r^(f,-f) 



gr(f,Q)-^t'""(f,o) . gr(f,o)-^r"(f,Q) 

= = ^ ^ 

g?'"(f,Q)-^t'""(f,o) , .0?'^(f,o)-c?'-^(f,o) 

= ■ = ± z = ^ 



— I 



x/-)l,a;/7r _ tt \ _/'/-)!, -x / tt _ tt \ /'/-)!, 2/ / tt _ tt \ _/'/-)!, -J/ / tt 



, ■ I4' 2/ '-^t ^.4' 2/ I I4' 2/ '-^t 



o^-(f,o)-o^-(f,o) . 0^^(f,o)-0^-^(f,o) 

P43(0 = 1 ±^ 1 

0?'-(f,Q)-0?'-(f,Q) . 0?'^(f,Q)-0?'-^(f,Q) 



■ IZ' 2^ W' 2/ ^ ^^4' 2^ \ai 



\l, — x 



+ 1 



— I 



)^,-yi 



2) ^4' 2^ _|_ '-^t I4' 2^ I 4 



4 ' 
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APPENDIX B 

After consistent absorption in the exponentials of terms linear in t, the entries Gij{t), 
i,j = 1, 2, 3, of the matrix Qt solution of equation (4.5), calculated up to first order in the 
dissipative term D, can be expressed as 



Quit) = e-^«+«)*cosa;ot + sinuot , {Bl) 

Qizit) = - — sin "^t (c cos "^t - Psm "^t) ; (53) 

Q2i{t) = fe-(«+«)* - — ) sincuot , (54) 

Q22 {t) = e- ("+«^^ cos uot + ^—^ sin uot , (55) 

LOq 

G23{t) = - — sin (p cos '^t + csin '^t] ; (56) 
UJo Z \ 2 2 / 



Gsi {t) = --sin^t(ccos^t + P sin ^t) (57) 
UJQ 2 \ 2 2 / 

032{t) = - Sin (c sin - Pcos ^t) (58) 

LOq 2 \ 2 2 / 

6?33(t) = e-2T* . (59) 
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